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O ! 

, We construct generally applicable small-loss rate expansions for the density operator of an open system. Suc- 
?H \ cessive terms of those expansions yield characteristic loss rates for dissipation processes. Three applications 
are presented in order to give further insight into the context of those expansions. The first application, of a 
two-level atom coupling to a bosonic environment, shows the procedure and the advantage of the expansion, 
t—{ t whereas the second application that consists of a single mode field in a cavity with linewidth k due to partial 
transmission through one mirror illustrates a practical use of those expansions in quantum measurements, 
' and the third one, for an atom coupled to modes of a lossy cavity shows the another use of the perturbative 
_ expansion. 

O ' PACS numbers:03.65.-w,05.30.Fk,42.50.Dv 

o . 

The study of open quantum systems has recently attracted the attention of physicists from various fields: 
cosmology[l], condensed matter[2], quantum optics[3-7], particle physics[8], quantum measurement [9, 10], and 

<*S i| quantum computational, 12]. The problem can be described generally as interest in the effective dynamics 
of one subsystem of several interacting subsystems. A formal framework to describe the effective dynamics 
of such a subsystem is set up in ref.[13], and a short-time perturbative expansion for coherence loss has also 

j— i been constructed[14]. To some extent(for example, if we are interested in a behavior for finite time), however, 
time is not as good as the loss rate as a perturbative parameter. Motivated by this and recent experimental 
developments [15-20] as well as the analysis of models related to them[21-24], we construct generally small- 
loss rate expansions for dissipation. The results suggest that these are useful in many areas such as high-Q 
Cavity QED[15-20], quantum computation [11, 12, 21-26], quantum measurement [27, 28], quantum optics[3-7] 
etc. 

We consider an open quantum system, the total Hamiltonian describing such a system is expressed as 



H — Ho + Henv + Hi, (1) 

where Hq and H env indicate the free Hamiltonian of the system and of the environment, respectively. Hj 
is the interaction Hamiltonian between the system and the environment. It is well known that the form of 
the master equation depends on the precise kind of the system-environment interaction. In order to derive 
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a master equation for a quite general Hi, let us suppose that, in the Schrodingcr picture, Hi can be written 

as [3,4] 

Hi = hJ2(X+A m + X-Al) (2) 

m 

where the X^axe eigenoperators of the system satisfying 

[H ,X±\ = ±hu> m Xt. (3) 

This form is quite general, since any system operator can be decomposed into eigenoperators of H - As 
shown in ref.[3,4], we can write the master equation in the following form (in the Schrodinger picture) 

p(t) = -i[H , P] + \Y, K m (2X mP X+ X+X-p P X+X m ) 

m 

+ GmC2X+pX- X-X+p P X m X+), (4) 



where 



POO 

K m = 2Re[ / dre^ T Ti env {A m (T)Al(0)p env }}, 
Jo 

r-OG 

G m = 2Re[ / dTe^ T Tr env {Al{T)A m {0)p env }} 7 
Jo 



p{t) = p(t, K m , G m ) stands for the density operator of the system and p env denotes the density operator of 
the environment. Notice from eq.(4) that G m should vanish at zero temperature T = 0, while K m should 
not if A m are indeed destruction operators of some kind. In case the constant G m and K m are smaller than 
any one of the internal coupling parameters of the system, the density operator may be expanded in powers 
oiK m and G m [2% 

p(t, K m , G m ) = p(t, 0, 0) + V ^-K m + V ^-G m 

m m 

1 x ^ d 2 p v 1^ d 2 p 

^ TZti 9K m dK n n 2 ^ dG m dG n 

m,n 

Substituting this expression into the master equation, we find the following set of equations 



p(t,0,0) = -i[H ,p(t,0,0)] (6) 
<h> = -i[ ff o, 4ft~\ + \i^+p{t, 0, 0)X m - X-X+p(t, 0, 0) - p(t, 0, 0)X m X+) 



dG m dG m 2 

9 P _ ,rr7 d P 1 , X /ov n ,,v- vH „ .„ „ mv 



= «^-] + n(2X-p(t, 0, 0)X+ - X+X mP (i, 0, 0) - p(t, 0, 0)X+X m ) (7) 



9X m '9if m J 2 

Q2 P - i\H dp 1 i l (1X- dp X+ X+X- dp dp X+X-) 



9x m 9x„ u, dK m dK n ' 2 V m aif„ m m m 9x„ dK n ' 
+ 1 (ix- dp x+ x+x- dp dp x+x-) 
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92 P = .rrr 9p ■ 1 + + dp ^Ly-V^ 

dG m 8G n [ °' dG m dGj 2 { m dG n m m m 8G n 8G n m m> 

+ \ {2X "^L X " X - X ™^L - ^L X - X " ] 

° 2 P j\H dp 1 I 1 (1X- 9P X + X+X- 9P 9P X+X-\ 



aif m ac„ u, dK m dG n > ' 2" m dG n m m m dG n dc n 

+ \ {2X "mL x " x ™ x ™^h ^t x - x - ] (8) 

Generally speaking, given an initial condition, p(0,0,0), we can solve eq.(6) exactly, which gives the zeroth 
order solution for the density operator p(i, 0, 0). Substituting the zeroth order solution into eq.(7), dp/dK m 
or dp/dG m can be calculated. Following this procedure, successive terms of the expansion (5) could be 
worked out, though the calculation is complicated. Some words of caution are now in order. From the 
mathematical point of view, the expansion (5) holds if and only if the series converges. This may be satisfied 
easily in physics for a large number of open systems. For example in a high-Q cavity, the loss rate of the 
atom-cavity system is small enough to permit us to expand the density operator in powers of the loss rate. 
At the end of this section, we will present some discussion about this point in detail. 

To illustrate the advantage of the expansion (5), we present here a simplest model, which describes a 
two-level atom coupling to a bose-mode environment. The master equation of such a system is given by[3] 

p = ~ifl[a z ,p] + \^j{2a^pa + - p<T + (T~ - cr + a^p} (9) 

with 

/>oo 

r Tr e „4MT)C(0)p e ™}], 



2wRe[ [ 
Jo 



oo 

a lU) m Tr 



'0 

where b^bm) stands for the creation(annihilation) operator of the m-th mode of the environment, Q is the 
Rabi frequency , and a z {u + , cr~) denote the Pauli matrices. To obtain the form of the master equation given 
in eq.(9), the environment was assumed to be in its vacuum state. According to eq.(5), (a z (t)} reads: 



- d P s , 1 2rp ( d 2 P 



(a z (t)} = Tr(p(i,0, 0)a z ) + jTr(^a z ) + - 7 ^Tr(— ^) + ... (10) 



The first term in eq.(10) is (a z (0)). In order to calculate Tr(|^a z ), we first evaluate Tr(§^a z ), it is given 
that 

Tr(|^) = -MO)). (11) 



Using the same procedure as mentioned, we arrive at 
The eqs.(ll,12) together give 



Tr ta^) = 2(^(0))*. (12) 



(a z (t)) = M0)) - Jt(a z (0)) + ^t 2 M0)) + ... (13) 
Using the algebra of Pauli matrices, we obtain straightforwardly from the master equation that 

Mt)> = M0)>e-*. (14) 
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A comparison between eq.(14) and (13) shows that for small 7, the expansions are a quite good approach 
for the two level dissipative system, and this result is quite general. 

Noticing the eq.(13) is expanded in the product of loss rate and time, we present here the other example 
to show that these expansions are generally in powers of the loss rate, but not in the product of time and 
the loss rate. Consider the master equation given in eq.(9) As mentioned above, we calculate (a z (t)} in order 
to illustrate the advantage of the expansions. The results of (a z (t)) show no difference between the short- 
time expansions and the small loss-rate expansions. To show the difference between the two expansions, we 
calculate (a x (t)}. For simplicity, we only present the results up to first order of 7. It follows from eq.(5) that 

(a x (t)) =Tr(p(t,0,OK)+ 7 TY(^a x ) + ... (15) 

It is easy to show that (setting U= 1, (a x (0)) = 1)) 

Tr(p(i, 0, 0)ct x ) = cos(m). 

In order to compute Tr(^a x ), we have to calculate Tr(^a x ). Based on the expansions, we arrive at 

Tr(|^) = Tr(-i[H Q , ^-<j x \) + \{2a~ p{t, 0, 0)a+a x - p(t, 0, 0)a+a~a x - a + a~p{t, 0, 0)a x ) = a + b. 



Simple calculation gives 



a = "" Tr( ^^ ) - 



As state above, in order to compute a, we have to calculate Ti(^a y ). Using the same procedure as above, 
we show that 

a = -n 2 J dtTr(^a x ) - (cos(ftf) - 1), 



and 



The results for a and b together give 



b = — cos(Oi). 



8 2 v(t\ 

' y[ ' +n 2 y(t) =2f2sin(ta), 



dt 2 

with y(t) = Tr(^a x ) and initial conditions y(t — 0) = 0,y(t)\t=o = —1- This is a two order differential 
equation and can be solved easily, once y(t) is known, (<r x (t)} up to first order of 7 is given. It is obvious 
that results given above are indeed different from the short time expansions, since the results given by short 
time expansions are in powers of time t. 

The results up to first order of 7(15) and an exact numerical results are illustrated in Fig.l. The parameters 
chosen are fl = 2, and time is in units of fi. In Fig.l the scattering line represents the exact numerical results, 
whereas the dot line and the solid line show the results from the expansion. The dot line and the solid line 
arc for different 7, and 7 for dot line is smaller than one in solid line, those curves show that the expansions 
(15) are indeed a good approximation to the exact solution. 
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We need to point out that for most open systems, average values of meaningful quantities can't be obtained 
exactly in any way. Therefore the expansions of the density operator provide a practical approach to the 
exact solution. For example, consider a single-mode field in a lossy cavity. The density operator for that 
mode obeys the following master equation in the Schrodinger picture[3,4], 

p = —i[ujfa^a, p] + —{2apa) — a^ap — pa^a), (16) 

where k is the linewidth of the cavity mode with frequency uif. In most textbooks, the solution of the 
master equation is given in terms of diagonal matrix elements (n\p\n) in a stationary state. Given an initial 
condition for the density operator, the evolution of p, however, is more useful than the stationary solution. 
In what follows, we present a solution of the master equation in a number state (Fock state) basis. 

For a high-Q cavity, the linewidth k due to partial transmission through one mirror is so small that we 
can expand p in powers of n: 

AO + J2.+ !££,■• + ..., (it) 



In a number state basis {\n),n = 0, 1, 2, 3...}, the expansion can be written as 

1 d k p m 



m.n m,n k—1 

Here the subscripts on the density operator p mn indicate matrix elements of p in the number basis and po 
is the solution of eq.(16) with k = 0. With this notation, it follows from eqs.(6),(7) and (8) that 

P°mn(t) = plMe- 1 ^- 11 ^, 

= Cfc (i) e -*"/(»- m )*. (19) 

This iterative equation gives the density operator expansions of the system under consideration. Here, p m „(0) 
stands for the initial condition of p, and 

Cfc(i) = f ' F k (t')e^ n - m ^' 'dt , 
Jo 

F 1 (t) = y/(n + l)(m+l)p° m+hn+1 (t) - \mp° mn {t) - ^np° mn (t), 

F k (t) = y/(n + l)(m + 1) Qj^> k = 2 ' 3 ' 4 - ( 2 °) 

The master equation in the form (16) is widely used in field-quadrature measurement[27,28]. As shown in 
Ref. [27,28], different approximations to eq.(16) correspond to different measurement schemes, therefore the 
expansions (6-8) for the density operator provide a new method to develop quantum measurement theory. In 
contrast to the short time perturbative expansions [14], the expansions (6-8) hold for finite time as long as 
the linewidth k is small. In other words, whether the expansions hold does not depend on time t. 

In addition to quantum measurement, these expansions have use in high-Q cavity QED. There are many 
interesting features in cavity QED. One of them is spontaneous emission. Spontaneous emission is so fun- 
damental that it is usually regarded as an inherent property of matter. The master equation for a single 
atom coupling to a mode of a lossy cavity is given in the interaction picture under rotating- wave and dipole 
approximations by [15] 
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Pi = — (2(7 pi<7 + — a + a pi — pja + <7 ) + —{2ap I a { — a) api — pja^ a), (21) 

where pi stands for the reduced density operator of the system that consists of an atom and a cavity mode, 
7 denotes the linewidth of the atom, and n describes the loss rate of the cavity. This is different from eq.(9) 
in which the loss of the single-mode field is neglected. It is well known that the emission spectrum may be 
expressed in terms of average values of the atom operator a. In this sense, we may calculate the average value 
of the atom operator to replace computing the emission spectrum without any loss of generality. Moreover, 
the study of many other effects in cavity QED such as atomic dipole squeezing [30], population trapping[31], 
and atomic collapse-and-revival phenomenon[32] may be reduced to calculate and analyse the average value 
of the atom operator. In the remainder of this paper, based on the expansion scheme, we compute the 
average value of an atom operator given by A = A( + V + + A( _ )<7 _ + \^a z . For this end, we first of all list 
the expansions of the density operator pi{t) in the interaction picture 

p l{t )=p I (t) + -g- K +— 1 +--^ r K +2^^ +2^T 7K+ - (22) 

Here, 

P°i(t) = 0, 

^ = \{2o- p»(t)a+ *+a-A{t) p?(t)o- + 0. 

^ = \{2ap°{t)cJ <Jap%t) p%tWa), 
d 2 Pi(t) I dpi + + _dp! dpi + _ 

-^r = 2 i2(T a G a } ' 

d 2 Pi(t) 1. dpi f f dpi dpi f 
-^=2 {2a ^ a - flB fe-& aB) ' 

d 2 Pi(t) 1 _dpi + + _dpi dpi + _ I dpi f f dpi dpi f 
^T=2 {2a to* a a ~d^ ~ ~d^ )+ 2 (2a a7 a a a 97-97° a) - 

It is easy to show that for any atom operator A, 

Tr(^A) = Tr(/^ A) = (23) 



for n > m ^ 0, while 



where _B = — 4\^a + — 4A(~)<j~. Eqs.(23) and (24) suggest that the average value for any atom operator 
can be calculated analytically as an expansion in powers of k and 7, provided p (t) (the zeroth order density 



6 



operator in the Schrodinger picture) is known. Generally speaking, given a initial condition for p, the p°(t) 
that obeys the von Neumann equation can be given readily. In the model presented above, the von Neumann 
equation is given by 

p° = -i[H ,p }, (25) 

Where H denotes the free Hamiltonian for the cavity-atom system (Jaynes-Cummings model) 

H a = ujja^a + -0J a a z + g(c^a~ + a + a). (26) 

If the cavity-atom system is initially in a state |e, n) — |e) ® i.e. the atom is in its excited state, while 
the single-mode cavity is in the number state |n), then p°(t) reads 

P (t) = \iP (t))(iP°(t)l (27) 

where 

|V°(*)> = isin^ +1 ( e -^+(" +1 )*- e -^-(" +1 )*)| 5 ,n+l) 

+ ( Sin 2 ^±l e -iE + (n+l)t + cos 2 ^±l e -iE.(n+l) t) ^ n y (2g) 

E±(n + 1) = ^(2n + 1) ± y/S 2 + g 2 (n + 1), 6 n+1 = arctg 2g ^~ ] ~ ,5 = u f -u a . (29) 
It follows from eq. (6,7,8) that 

Tr(p°(t)A) = ±Tr(p Q (t)B) + A<*> ]T{| sin 2 ^I e -«M»+i)* + cog 2 ^i e -iB_(n+i)t|2 

n 

- sin 2 6 n+1 sin 2 (^S 2 +g 2 (n + l)t)} 
Tr{Bp°(t)) = ^{-2A« sin0„(e- iB +(") t - e -^-(™)*) 

n 

( sin 2 e _^+L e iE + (n+l)t + CQS 2 iE-( n +l)ts _ 

v 2 2 ; 

2A(-)(sin 2 ^+L e -iE + (n+l)t + CQS 2 «n+l - iB _(n+l) t) 

v 2 2 ' 

sin^^e^+W'-e^-W')}- (30) 

Then successive perturbative terms of average value up to second order of 7 for an atom operator are given 

by 

(A)(t) = Tr(p (t)A) + lj\(p°(t)B)dt -f-J*dt'J* Tr(p°(t")B)dt" + .... (31) 

Based on the short-time expansion, (A)(t) = Tr(Ap(t)) ~ A + Ait + A 2 t 2 + in powers of t, where 
A Q = Tr(Ap(0)), A 1 = Tr(A§f (0)) and A 2 = ±Tr(,4§^(0)). This is quite different from the results given by 
eq.(31). 

Although we are currently investigating the perturbative expansion for an open system, we opt here for 
a few qualitative comments. Mathematically, the perturbative expansion is a good approach to the exact 
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solution of the master equation so long as loss rates 7 and k are smaller than all other internal coupling 
constants of the system. This condition holds for high-Q cavities from the physical point of view. In fact, 
an optical cavity of ~ 20pm diameter has g/2n ~ \2hMHz and k/2it <~ lOOKHz for reasonable Q ~ 10 9 . 
Thus the ratio g/n ~ 10 3 . Even g/n ~ 10 4 seems feasible for microspheres [33]. Generally, in the optical 
domain 3/7 ~ 10 2 , great enough for the perturbative expansion in powers of 7 to hold. 

In the end of this paper, we turn our attention to study the decoherence in N two-level atoms. This 
problem is usually related to the register in quantum computer. A few papers[ll,12,34] have been published 
on this subject , but a key additional feature of the present paper is to study the decoherence from a new 
aspect. If the system consists of N two-level atoms, the decoherence is due to the inevitable coupling of the 
N atoms to the external environment. Generally, the environment may be treated as that consists of an 
infinite number of oscillators. The Hamiltonian describing such decoherence process takes the form 

H = H s + H env + Hi , 

N 

H s = ^ Ojof, 

i=l 

oc 

H env — ^kbfrbk, 
fc=i 

N 00 

H I = Y,Y,^ub\c7i +h.c), (32) 

i=l k=l 

where erf are the spin-i Pauli operators (i denotes the qubit index) and bk,b\. are the bosonic operators, 
H s , H env are the free Hamiltonian of the system and the environment, respectively. And Hi stands for 
the N qubits-environment interaction. This model is closely related to the Dicke maser model[35,36]. The 
Hamiltonian (32) is complicated so it is hard to find its exact solution though the Hilbert space associated 
with this model can split into invariant eigenspaces[34]. Fortunately, with the perturbative approach created 
in the previous section, the complex system can be easily treated. To start with, we give the master equation 
of the system 

p = -i[H, p] + Ki{2arpaf - pafar - a+a^p)} 

i 

+ Gi(2afpar - a^afp - P <yr<j+)} 

i 

= -i[H,p]+Cp. (33) 



Here 



K m = 2Re J dTe in ™ T TY env {A m (T)Al(Q)p env }, 
G m = 2Re J dTe in ™ T Tr env {Al(T)A m (0)p env }, 



A m (r) = ^9m j b j e-^ T . 
i=i 
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In the case of so-called Dicke limit [35,36], A m does not depend on the atom index to. This holds, for example, 
when the typical environment wavelengths are much greater than the distances between the atoms. 
In order to study the decoherence of the atoms, we assume that the initial state of the system is 

N 

|VU=S+|0)+ [] 1°))' ( 34 ) 

j=m+l 

where S+ = UT=i N cr/|0), and |0) = |0)i ® |0) 2 ® ... O |0}jv stands for the lower state of the atoms. Eq.(34) 
indicates that there are m atoms in the upper state |1), and the rest of the N atoms are in their lower state. 
With those initial conditions, the probability of the atoms remaining in the initial state is given by F m (t) 

F m (t) = l--^-^--^--^-^ + 0( G 3) + 0(K% (35) 

1 K 1 J- if 2 1 Gl 1 G2 i-GK 

with 

1 



= t[2 + 2(2m- N)], 



Kl 



^— = —t 2 [2 + 2(2m — N)], 

K2 



-L=t[2-2(2m-N)}, 

1 Gl 



-}- = -t 2 [2-2(2m-N)], 

1 G2 



1 -GKt 2 [A(2m- N) -4]. 



T GK 2 

Here, we suppose that all qubits are alike, so TO = £1 and K m = K and G m — G. G and K = G + 1 depend 
on environment temperature T through G = l/[exp(KQ/kT) — 1], which indicate that the probability decrease 
with the temperature increasing. In fact, the fidelity in the field of quantum information is nothing but an 
overlap between the initial and final state of the qubits(two-level system). Eq. (35) suggests that the fidelity 
depends on to i.e. the number of the atoms in upper state initially. And to get the maximum of the fidelity, 
the variable to should be taken as small as possible. 

To sum up, in this paper, we construct the small-loss rate perturbative expansion for the density operator 
of an open system. The expansions provide a quite good approach to the exact solution in case the master 
equation of the system can not be solved exactly. As an interesting application of this expansions, we used 
it to calculate some average values such as a z and a x in dissipative two-level system, the expansions of the 
density operator for a single-mode field in a lossy cavity are also presented, and the dynamical property in 
N two-level atom system. 

In addition, the other meaningful quantities of the open system such as energy, occupation probability etc. 
can be expanded in the same spirit of the density operator, so long as the master equation of the system is 
known. 



9 



[1] M.Gell-mann, J.B.Hartle. Phys. Rev. D 47(1993)3345. 

[2] A.O.Caldcira and A.J.Leggett, Ann. Phys.(N.Y.)149(1983)374. 

[3] C.W.Gardiner, Quantum noise (Springer- Verlag), New York, 1991. 

[4] William H.Louisell, Quantum statistical properties of radiation, John Wiley, 1973. 

[5] Xue-Xi Yi,Physica scripta,52(1995)508. 

[6] Xue-Xi Yi, J.Mod. Opt. 43(1996)675. 

[7] Xue-Xi Yi, J.C.Su J. Phys. B 30(1997)L411. And references therein. 

[8] H.T.Elze, Nucl. Phys. B 436(1995)213. 

[9] W.H.Zurek, Phys. Rev. D 26(1982)1862. 
[10] D.F.Walls, M.J.Collet, G.J.Milburn, Phys. Rev. D 32(1985)3208. 
[11] W.G.Unruh, Phys. Rev. A 51(1995)992. 

[12] R. Landauer, Philos. Trans. R. Soc. London Ser. A 353(1995)367. 
[13] M.C. Nemes and A.F.R.de Toledo Piza, Physica 137A(1986)367. 

[14] Ji II Kim, M.C. Nemes, A.F.R. de Toledo Piza, H.E. Borges, Phys. Rev. Lett. 77(1996)207. 
[15] W.Ren and H.J.Carmichael Phys. Rev. A 51(1995)752. 
[16] S.Haroche and D.Kleppner, Physics today 42(1989)24. 
[17] P.Meystre, Phys. Rep. 219(1992)243. 

[18] O.Benson, M.Weidinger, G.Raithel, and H.Walther, J.Mod. Opt. 44(1997)2011. 

[19] M.S.Kim, G.Antesberger, C.T.Bodendorf, and H.Walther, Phys. Rev. A 58(1998) R65. 

[20] D.W.Vernooy, A.Furusawa, N.Ph.Georgiades,V.S.Ilchenko, H.J.Kimble, Phys. Rev. A 57(1998)R2293. 

[21] J.I.Cirac, P.Zoller, Phys. Rev. Lett. 74(1995) 4049. 

[22] S.J.Van Enk, J.I.Cirac, and P.Zoller, Phys. Rev. Lett. 78(1997)4293. 

[23] J.I.Cirac, P.Zoller, H.J.Kimble, and H.Mabuchi, Phys. Rev.Lett. 78(1997)3221. 

[24] S.J.Van Enk, J.I.Cirac, P.Zoller Phys. Rev. Lett. 79(1997)5178. 

[25] A special issue on coherence and decoherence, Proc. R. Soc Lond. A 454(1998)257(begin page). 
[26] C.H.Bennett, Physics Today, 48(1995)24. 

[27] V.B. Braginsky, F. Y. Khalili, Quantum measurement (Cambridge, 1992). 
[28] H.M.Wiseman, G.J.Milburn, Phys. Rev. A 47(1993)642. 
[29] N.G.van Kmpen, Phys. Rep. 124(1985)69. 
[30] T.Nasreen, J. Opt. Soc. Am. Bll(1994)386. 

[31] B.D.Agap'ev, M.B.Gornyi, and B.G.Matisov, Yu.V.Rozhdestvenskii, Phys. Uspekhi 36(1993)763. 
[32] P.R.P.Jushi, S.V.Lawande, Phys. Rev. A 44(1991)2135. 

[33] M.L.Gorodetsdy, A.A.Savchenkov, and B.S.Ilchenke, Opt. Lett. 21(1995)453. 
[34] P.Zanardi and M.Rasetti, Phys. Rev. Lett. 79(1997)3306. 
[35] R.Davidson, J.J.Kozak, J.Math. Phys. 8(1971)903. 
[36] K.Hepp, E.Lieb, Phys. Rev. A8(1973)2517. 

Figure captions 

Fig.l: a x (t) vs time t, the parameter chosen is O = 2. The scattering line represent the exact numerical 
results, whereas the solid line and the dot line show the results from the expansion, 7 in the solid line and 
the dot line is different, dot line:7 = 0.01, solid line:7 = 0.05. 
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